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<^ Abstract 

A reduction of the Maxwell's system to a Fredholm second-kind integral equation with 
weakly singular kernel is given for electromagnetic (EM) wave scattering by one and 
many small bodies. This equation is solved asymptotically as the characteristic size 
of the bodies tends to zero. The technique developed is used for solving the many- 
body EM wave scattering problem by rigorously reducing it to solving linear algebraic 
systems, completely bypassing the usage of integral equations. An equation is derived 
for the effective field in the medium, in which many small particles are embedded. A 
method for creating a desired refraction coefficient is outlined. 
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jj^ '. 1 Introduction 

ON 

There is a large literature on electromagnetic wave (EM) scattering (pQ- [3]), to name a few 

books. In [5] (see also [6], [7] and references therein), wave scattering theory is developed 
for small bodies of arbitrary shapes. The ideas and methods from [5] were used in recent 
papers [8], [11], where a many-body scattering problem was solved asymptotically for 
scalar wave scattering, as the characteristic size a of small bodies tends to zero. One of 
the applications of this theory in acoustics is a method for creating new materials with the 
desired refraction coefficients. In particular, the new materials may have a desired wave- 
focusing property: a plane wave incident at a fixed direction with a fixed wavenumber will 
be scattered by the new material with a desired radiation pattern. Another application 
([13], (p3|)) is a method for creating materials with negative refraction, that is, materials 
in which the group velocity is directed opposite to the phase velocity. 

In [13] the new materials are created by embedding small particles into a given material 
in a bounded domain D so that the number of the particles embedded into a small subcube 
A C D, is given by the formula 

Af(A) = -^ / N(x)dx[l + o(l)], a -+ 0, (1) 

a J A 
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where a is the radius of the embedded spherical particles D m = B(x m , a), x m is the center 
of the ball B(x m ,a), a is its radius, < x < 1 is a parameter, N(x) > is an arbitrary 
continuous in D function, which we can choose as we want, the boundary condition on 
the surface S m of the m-th particle D m is: 



du h(x m ) 

dN = ^^ U ° n5m ' (2) 



where N is the unit normal to S m , pointing out of D m , and h(x) is an arbitrary continuous 
function which we can choose as we want. 

The purpose of this paper is to develop a theory for the electromagnetic (EM) wave 
scattering, similar to the one developed in [8]- [12], [14], [15], in the sense that it yields a 
method for creating materials with a desired refraction coefficient. Condition ([2]) is not 
used in this paper, but formula (f73j) of this paper is a generalization of ([T|). 

In Section [2] we present a theory for EM wave scattering by a single small body. In 
Section [3] many-body scattering problem is solved. In Section [4] the negative refraction is 
discussed. In Section [5] the conclusions are formulated. 

2 EM wave scattering by a single small body 

Let D be the body and a := | diamD. Assume that 

ka < 1, (3) 

where k > is the wavenumber. The governing equations are: 

V x E = iujfiH, V x H = -iu>e'(x)E inM 3 , (4) 

where uj > is the frequency, \x = const is the magnetic constant, e'{x) = e > in 
D' = R 3 \ D, e'(x) = e(x) + i^-; a(x) > 0, e'(x) / Vx G M 3 , e'(x) G C 2 (M 3 ) is a twice 
continuously differentiable function, and cr(x) = in D' is the conductivity. From (4) one 
gets 

V x V x E = K 2 (x)E, H= — - — , (5) 

K 2 (x) := to 2 e'(x)fi. (6) 

We seek the solution of the equation 

VxVx£ = K 2 (x)E (7) 

satisfying the radiation condition: 

E(x) = E (x) + v, (8) 

where Eq(x) is the plane wave, 

E (x)=£e* ka - X , k = -^ (9) 



where c = uj-y/ejl is the wave velocity in the homogeneous medium outside D, e = const 
is the dielectric parameter in the outside region D', \x = const, a £ S 2 is the incident 
direction of the plane wave, £ ■ a, = 0, £ is a constant vector, and the scattered field v 
satisfies the radiation condition 

dv , ,1s 

— ikv = o(-), r = \x \ — > oo, (10) 

or v r y 

uniformly in directions (3 := f. 

If E is found, then the pair {E, H}, where H = , solves our scattering problem. 

Our goal is to derive a Predholm second-kind integral equation for E. 
Let us rewrite equation ([7]) as 

-AE + V(V-E)-k 2 E-p(x)E = 0, (11) 

where 

p(x) := K 2 (x) - k 2 , p(x)=0mD'. (12) 

Note that ImK 2 (x) = Imp(x) > 0. This will be used in the proof of Claim 3 below (see 



It follows from that 

= V • (K 2 (x)E) = VK 2 (x) • E + K 2 (x)V • E. (13) 

From (fT2|) and (fT3|) one gets 

- A^ - k 2 E - p(x)E - V{q{x) -E) = 0, (14) 

where 

. N VK 2 (x) . 

q(x) ■= R Jy q(x) = Oin D'. (15) 

From (fT4"|) and (JSJ) one gets 

E = E + J g(x, y) (p(y)E(y) + V y (q(y) ■ E(y))\ dy, 



(16) 



where 



& ik\x-y\ 



4tt\x — y\ 

If q vanishes at the boundary of D, i.e., X7K 2 (x) = X7p(x) vanishes at the boundary of D, 
then, after an integration by parts in the last term of (|16p . and taking into account that 
V x g(x,y) = —V y g(x,y), one gets an equivalent equation, which we also refer to as (fT6|) : 

E = E + / g(x,y)p(y)E(y)dy + V x / g(x,y)q(y) ■ E(y)dy. 
Jd Jd 



Since we have assumed K 2 (x) ^ and K 2 (x) £ C 2 (M 3 ), it follows that equation (|16p 
is a Fredholm equation of the second kind, because the integral operator T in (|16p is 
compact in H l {D), where H l (D) is the usual Sobolev space. 



Indeed, the operator Bf := f D g(x,y)f(y)dy acts from L 2 (D) into H 2 (D), and the 
operator B\E := p(y)E + V(g(y) ■ E) acts from H 1 ^) into L 2 (D). Thus, T acts from 
H l {D) into H 2 (D) and is, therefore, compact by the embedding theorem. 

Let us summarize the result. 

Lemma 1 // K 2 (x) G C 2 (R 3 ), K 2 (x) / 0, K 2 (x) = k 2 > in the exterior domain 
D' : = BL 3 \ D, then the operator T in (|16p is compact in if 1 (Z)) ; so that equation (|16p zs 
of Fredholm type in H^^D). 

Lemma 2 Equation (|16p is uniquely solvable in H l (D). 

Proof. Let us prove Lemma 2. It is sufficient to prove that the homogeneous version 
of equation (|16p has only the trivial solution. If E solves the homogeneous equation (|16p . 
then E solves equation (|7|) and satisfies the radiation condition (|10p . The only solution to 
([7]) satisfying (|10p is the trivial solution E = 0. 

Let us give details and prove the above claims. 

Claim 1 . A solution to (fT6|) satisfies equation (J7J), and, consequently, (fT3|) . that is, 
V • (K 2 (x)E) = 0. 

Thus, Claim 1 states that (fT6|) is equivalent to the original equation ([7]), which is not 
at all obvious. 

Claim 2 . The only solution to fl7J), (HUD is S = 0. 

Proof of Claim 1. If E solves (1161) and is of the form (TSl) , then 



(-A - k 2 )E = p(x)E + V(q(x) ■ E). 
Rewrite this equation as 

VxVxB- V(V • E) - V(q(x) ■ E) = K 2 (x)E, 

or 

„ „ ^ ^(K 2 (x)S7-E + VK 2 (x)-E\ „. Nr ^ 

VxVx£-V \—U R2{x) { ' j = K 2 (x)E (18) 

Denote 

-±-^V-(K 2 (x)E):=iP(x). 

Taking divergence of equation (fTHj) one gets 

- A^ - K 2 (x)^ = in M. 3 . (19) 

The function if) satisfies the radiation condition (1101) . 

Claim 3 . The only solution to (|19p . which satisfies condition (|10p . is ^ = 0. 

We prove Claim 3 below. Assuming that this claim is proved, we infer that ip = 0, so 
V • (K 2 {x)E) = 0, and equation © holds. 

Claim 1 is proved. □ 

Let us prove Claim 3. Equation (|19p can be written as: 

- Aip - k 2 tp - p(x)ip = inK 3 , (20) 



where p(x) is defined in equation (I12|) . and p(x) = in D'. It is known that if k 2 > 
0, p{x) 6 L 2 (D), Imp(x) > 0, D is a bounded domain, and ^ satisfies the radiation 
condition (fTOl) and equation (j20j) . then ■*/> = (see, e.g., [8]). 

For convenience of the reader, we sketch the proof. From (|20p and its complex conju- 
gate one derives the relation 

ipAip - ipAtp - 2ilmp(x)\ip\ 2 = 0. (21) 

Integrate (j21|) over a ball Br centered at the origin of radius R, and use the Green's 
formula to get 

/ U-$- -$-$-) ds-2i / Imp(x)|V'| 2 dx = 0. (22) 

J\x\=R dr dr> J D 

By the radiation condition (|10p for ^ one can rewrite (j22j) as 



2ifc lim 

fl-»oo./| a .| =fl j D 



|V>| 2 ds-2i / Imp(x)|V>| 2 cte = 0. (23) 

JD 

Since Imp(x) > by our assumptions, and k > 0, it follows that 

lim / \ifj\ 2 dx = 0. (24) 

R-*°°J\x\=R 

Relation (|24p and the equation 

(A + A; 2 )V' = |x|>,Ro, (25) 

where Br D D, imply ip = for |x| > i?o- (See [7, p. 25]). 

By the unique continuation principle for the solutions of the homogeneous Schrodinger 
equation (|20p . it folows that ip = 0. Claim 3 is proved. □ 

Let us prove Claim 2. This will complete the proof of Lemma EJ 

If E solves ([7]) and satisfies (fl~0|) . the the pair {E, H}, where H is defined by the 

second formula in ([5]), solves the homogeneous Maxwell system ([!]), satisfies the radiation 

condition, and {E, H} is in i?^ (R 3 ). It is known (see, e.g., [3]) that this implies E = 

H = 0. Lemma [2] is proved. □ 

From formula (|16p assuming that a := ^diamD is small, ka <C 1 and the origin is 
inside D, one gets: 



E = E Q + e —±- 

\X\ 47T 



as \x\ — * oo, A := /3, and 



/ e-^(p(y)£(y) + V(?(y) • E))dy[l + O(^)] (26) 

JD \ x \ 



e- ik(3 - y = l + 0(ka), \y\<a. (27) 

Thus, 

E(x) = E (x) + ^ TT — (p(y)E(y) + V(q(y) • £?))d|/[l + 0(^-)}. (28) 

p| 47r JD \ x \ 



Here we have used the estimate 

i \v\ 2 

(\x\ 2 - 2x ■ y + \y\ 2 ) 2 = \x\ - (3 ■ y + 0(-f+ r ), b|>a>|y|, /3:=x/\x\ 

\x\ 

which yields 



ik\x-y\ Jk\x\-ik0-y+O(-£ r ) Jk\x\ _ u n : 



4vr|x-y| 4vr|x| (l + O(jfj)) 4?r 

where |y| < a < |x|. If &a < 1, then O(tjx) <C 0(t§t), and we get estimate (J2SJ). Usually 
the error term in (|28p is written as O(m) as |x| — ► oo, but for small a the term O(A) is 
small already for \x\ ^> a, and x may be still small. 

If D does not contain the origin and x m £ D, then one can rewrite fll6|) as: 

S(x) = E (aO+#(x,x m ) / p(y)E(y)dy + V x g(x,x m ) / q(y) ■ E(y)dy 



+ / [g(x,y) - g(x,x m )]p(y)E(y)dy + V x [g(x,y) - g(x,x m )]q(y) ■ E(y)dy, 
Jd ' Jd 

(29) 

where g is defined in (|17|) . Let 

^ m := / p(y)E(y)dy, u m = f q(y) ■ E(y)dy, (30) 

where p{x) is defined in (|12|) and q(x) is defined in (|15p . 
If 

d := |x — x m | ^> a, a := - diamD, 

then 

k 1 
\g(x,y) - g(x,x m )\ = \V y g(x,y) ■ (y -x m )\ < eamax(-, -^), (31) 

where c > stand for various constants independent of a, y is an intermediate point, and 
d = \x — x m |. Thus, 



Jim. •" 



0(z,2/) - g(x,x m )]p(y)E(y)dy 

D 



k 1 
< camax(-,-2)|p(x m )E'(x m )||L>|, 



where \D\ = 0(a 3 ) is the volume of D, and we assume that p(x) and E(x) are continuous 
inD. 

On the other hand, 

0( ^ Xm)E{ d Xm)m ) >> J m , (32) 



g(x,x m ) I p(y)E(y)dy 

D 



provided that | -C 1 and ka -C 1. These assumptions hold throughout the paper 



We have checked that 
Similarly, 
where 



•Jm **^ -'m- 



J' m <C l' m , a <C d, ka <1, 



/' := 



V x g(x,x m ) / q(y) ■ E(y)dy 
D 



J' - 



Also, 



V x / [ff(x,y) -3(a;,a; m )]q(y) ■ E{y)dy 

D 



Jm ^ -/ m , J m ^ ±rn- 



(33) 

(34) 
(35) 



Therefore, in the region \x — x m \ 3> a one can calculate E{x) by the approximate 
formula 

E(x) = E (x) + g(x, x m )V m + S7 x g(x, x m )z/ m , (36) 

with the error 0(^ + fca). Here V^ and v m are defined in (f30j) . 
Let us derive formulas for the vector V m and the scalar v m . 
Multiply (|2"9"j) by p(x) and integrate over D to get: 



V m = V m+ / p(x)5(x,a; m )dxKi+ / p{x)V x g(x, x m )dx v m , (37) 

Jd Jd 



where 



Vom ■= / p{x)E (x)dx, 
Jd 

and we have neglected the two last terms in (|29p . These terms are smaller that the ones 
we kept, as we prove soon. 
Let 

/ p(x)g(x,x m )dx := a m , / p(x)V x g(x,x m )dx := A m . (38) 

Jd Jd 

Take dot product of (|29p with q(x), integrate over D and again neglect two last terms to 

get 

V m = ^Om + B m ■ V m + bmVm, (39) 

where 

(40) 



VQm ■= / q(x) ■ E (x)dx, 
Id 



B 



l H (x)g(x,x m )dx, b m := / q(x) -V x g(x,x m )dx. 



D 



From (|37j) one gets 



£> 



m = 71 ^ 7T_ Vn 

1 #m 1 Q>rn. 



(41) 
(42) 



Note that if the body D is small, then \a m \ < 1, and actually \a m \ <C 1, so (142p is well 
defined. Also, \b m \ < 1 if D is small. From ([42]) and (|39|) one obtains: 

(1 — o, m )uo m + B m ■ Vom . . 



(, J- ^m J V J- "m J i>m " -™-m 

and 

t^ _ »0m . -™m (.1 — a m) l/ 0in + -Dm ' »0m f/M"\ 

J- "in J- ^m V-l- ^m/v-*- "mj "m ' -"-m 

Formulas (f36j) . (|4"3"]l . and ([4"4"j) solve £/ie scattering problem of EM waves by one small body 
D located so that x m G D. 

Let us now check that one may neglect the terms we have neglected. 

For example, one has 



D 



dxp(x) [g(x,y) - g(x,x m )]p(y)E(y)dy = / dyp(y)E(y) dxp(x)[g(x,y) - g(x,x m )] 
d Jd Jd Jd 

dyp{y)E(y)[ip 1 (y) - ipi(x m )] 

< / dy\p(y)E(y)\, 
Jd 

because that ipi(y) ■= f D dxp(x)g(x,y) is a continuous function. 

The function ipi(y) is continuous if p(x) G L 2 (D), since t^i G H 2 oc (R. 3 ) an< ^ the embed- 
ding from H 2 (D) into C(D) is continuous (and compact if -D C M 3 is a bounded domain, 
which we have assumed). 

A similar argument holds for tp2{y) '■= f D [V x g(x, y)—V x (g{x, x m )] ■q(y)dy and justifies 
the possibility to neglect the second term we have neglected. 

To summarize: 

If p(x) and q(x) are continuous functions then the scattering of EM waves by a small 
body D is solved analytically by the formulas (|36p . (|43p . (|44p . with an error of order 
0(% + ka). 

3 EM wave scattering by many small bodies. Creating ma- 
terials with a desired radiation coefficient. 

The basic equation now has the form 

M . M 



E{x) = E Q (x) + ^2 9(x,y)p(y)E(y)dy + ^ V x / g{x,y)q{y) ■ E(y)dy 

m=1 JD m m=1 J D m 

M M 

= E (x) + ^2 [g( x >%m)V m + V x g(x,x m )u m ] + y^(J m + K m ), 



m=l m=l 

where M is the number of small bodies, 



J m -= [g(x,y) - g(x,x m )]p(y)E(y)dy, (46) 

JD m 



and 

K m :=V x [g(x,y)-g(x,x m )]q(y)-E(y)dy. (47) 

JD m 

We assume throughout, that the points x m € D m are distributed in D so that their 
number A/"(A) in any subdomain A G D of D is given by formula (|73p (see below) with 
<p(a) = a 3_>< , xr > (see Lemma 3 below). As a — > and M — > oo, the set of the points 
x TO becomes everywhere dense in D. 

Neglecting J m and K m , let us derive a linear algebraic system (las) for finding V m and 
v m . If V m and f m , 1 < m < M, are found, then the EM wave scattering problem for M 
small bodies is solved by the formula 

M 
E(x) = Eq(x) + ^2 [g(x, x m )V m + V x g(x, x m )u m ] (48) 

m=l 

with an error 0(% + fco) in the region mini< m <jv/ \x — x m \ := d ^> a. 

To derive a linear algebraic system for V m and v m , multiply ([43]) by p(x), integrate 
over Z)j, and neglect the terms J m and K m , to get: 

M 

^ = V 0j + Yl ( a i™ V ™ + B jm v m ) , 1 < j < M, (49) 

m=l 

where 



Vqj := / p(x)E (x)dx, a jm = / p(x)g(x,x m )dx, B jm := / p(x)V x g(x,x m )dx. 

JDj J Dj JDj 

(50) 
Take the dot product of ([45]) with g(x), integrate over Dj, and neglect K m , J m , to get: 



M 
V 3 = U 0j + /__, [CjmVm + dj m U m ) , (51) 

771=1 

where 

7y j = / q(x) ■ E (x)dx, C jm := / q(x)g(x,x m )dx, d jm := / q(x) -V x g(x,x m )dx. 

JDj JDj JDj 

(52) 
Equations (|49p and (|5ip form a linear algebraic system for finding V m and 7y m , 1 < m < M. 
This linear algebraic system is uniquely solvable if ka <C 1 and a <ti d. The elements Cj m 
and -Bj m are vectors, and aj m , dj m are scalars. Under the conditions 

M 

max ^2 (\a jm \ + |dj m | + ||Bj m || + ||Cjm||) < 1, (53) 

— — m=l 

one can solve linear algebraic system ()49|) and (|5ip by iterations. In (|53p ||-B,,- m || and ||C.,- m || 
are the lengths of the corresponding vectors. Condition (|53p holds if a <C 1 and M is not 
growing too fast as a — ► 0, not faster than 0(a^ 3 ). 



Consider the limit, as a — ► and M — ► oo, of the field i£(x) := -E(x, a). Let us assume 
that the number /V(A) of small particles in any open subset A C D is given by the 
formula: 

A/"(A) = -^- / iVOr)^ [1 + o(l)] , a -► (54) 



where < xi < 1 is some number, and iV(a;) > is a continuous in D function. If ([57 
holds, then the distance between two neighboring small particles is of the order 0{a yci ). 
We choose k\ = (3 — x)/3 above formula ([73]) . 
Let us write (H31) as 



#0*0 = E Q {x) + S~] g(x,x m ) p(y)dyE(x m ) + V x g(x,x m ) y ■ E(y)dy 

(55) 
where the terms J m and K m are neglected. 

Choose the function p so that suppp = U m=1 2? m , where Z? m := B m := {y : \y — x m \ < 
a}, that is, the small particles are balls of radius a centered at the points x m , and suppp 
is the support of p. Choose p of the form: 

p(y)=p(r,a), y £ B m , r := \y - x m \, 1 < m < M, (56) 

where the function p(r, a) should be chosen so that 

/ p(y)dy = 0(a s ), a —* 0, s := 3x\ = 3 — x, (57) 

where < x < 3. One may like to have Z m = 0(a s ) as a — > 0, where 

with the same s as in (]58|) . If this would be possible, then it would make possible passing 
to the limit as a — > similarly to the passing done in p3]. However, as we will see below, 
the relations (]57p and Z m = 0(a s ) are not compatible, in general. 
One has: 

w I \ If \ V ~ x m. I dp . . 

Vp(y) =p(r,a)r u , r u := , p := — . (59) 

r ar 

Here r° is a unit vector and r is the length of the vector y — x m . Let 

^(y) = E{x m ) + rej — |- — r, u + 0(a ), r( := , (60) 



Cartesian basis of M 3 , and Sij 



(x m )i := x m ■ a, a ■ b is the dot product of vectors a and b, ei ■ ej = 5« is the standard 

0, * + j, 

1, i = 3 ' 
Over the repeated indices in (]60p and everywhere below summation is understood. 



10 



If p(y) = p(r, a) then 

w £(*„,) 

B m k 2 +p(y) 

for any constant vector E(x m ) by symmetry. 
Thus 

By symmetry the integral in (1621) vanishes if i ^ j. This integral for i = j is equal to 

rd^rpp^ r d0s[n9 j 2 * d 2 i = 123 r o. (63) 

J k 2 +p(r,a) J J 

One has r| = cos 2 0, r| = sin #sin 99, r 2 = sin 6 cos 2 (p, so 



/■27T /j 

<i6> sin 6> / d^r^ = — , 
Jo 3 



and, with r = to, one gets: 



4tt a f 1 dtt 3 p'.(at,a) , 4 . ., „ „ ,„,. 

Yi = —a A k 2+ \ a ta) +°( a )' a ^ ' * = 1 > 2 ,3. (64) 



Thus 



4tt j w ,„ „, , , . /" TjL(at,a)t dt 
Z m ~ Y a*I(a)V-E( Xm ), m .^ J l£-U— r = ta, (65) 



and 



im := / pdy = 47ra / dtt p(at,a) := 47ra j(a). 

JB m Jo 



Let 



7m 



(1 - t) 2 h(t), < t < 1, t := §, x = const > 0, 



a 3 Z" 1 
p(y)cfy = — 7m / (1 - t) 2 h{t)t 2 dt := c ln 

a Jo 



(66) 



P(r)=p(r,a) = i 4 ™^ y w ' ^ ^ > ^ a ' ' (67) 

where 7 m , Im7 m > 0, is a number we can choose at will. Then 



j m ■■= / P(y)dy = —7m / (1 - t) 2 h{t)t 2 dt := c lm a 6 - >t , (68) 



where c\ m is proportional to 7 m , so we can also choose e im at will. By (|65p . 

I{a)=l ltZ{lm{ \: t? T'^ v-=^- (69) 

7o ^ + 7m(l-0 2 ^(i) 

Can one choose /i(t) so that 1(a) = Ca~ ie [l + o(l)] as a — ► 0? If yes, then Z m and j m are 
of the same order as a — ► 0. 

11 



In fact, it is not possible to choose such an h(t). One has 

I' 1 

I(a) = t 3 ln[v + lm (l-t) 2 h(t)]\l-3 t 2 \n[v + lm (l-tfh{t)]dt. (70) 

J u 

Thus 

l(a)=lm/ + 0(l) = xlna + 0(l), a -> 0. (71) 

Therefore no choice of h € C[0, 1] can ensure the asymptotics 1(a) = Ca~ >c [l + o(l)] 
as a — ► 0. 

We have obtained the asymptotic formulas: 

jm = ci m a 3 -", Z m = ^(V-E)(x m )a 3 lna[l + o(l)], a — 0. (72) 

Let us take /i(i) = 1. Then c\ m = 7 m /30. To pass to the limit a — ► in formula ([55]) we 
use the following lemma, which a generalization of a lemma from [8] : 

Lemma 3 If f € C(D) and x m are distributed so that 

MA) = — l-r /" N(x)dx [1 + 0(1)1, a -*•(), (73) 

V(a) Ja 

/or any subdomain Acfl, where (f(a) > is a continuous, monotone, strictly growing 
function, ip(0) = 0, i/ien 

lim V/(x m )^(a) = / f(x)N(x)dx. (74) 

Proof. Let D = U p A p be a partition of Z? into a union of small cubes A p , having no 
common interior points. Let |AJ denote the volume of A p , 5 := max p diam A p , and y be 
the center of the cube A p . One has 

SE^ m ^ = S Y /(y (p) )]™^( a ) Yl l 



hm V/(y^)iV(yW)|A p | = / f(x)N(x)dx. 



The last equality holds since the preceding sum is a Riemannian sum for the continuous 
function f(x)N(x) in the bounded domain D. Thus, Lemma 3 is proved. □ 

Let <p(a) = a 3 '*, x > 0. Then, by formulas ([55]). ([65]) and (|72|). one gets: 

lim —P~ = 0, lim — ^- = c im / 0. (75) 

a— >0 (f{a) a^0 ip{a) 

Therefore, by Lemma [3] the right side of formula (j55f) has a limit as a — ► provided 
that (|54j) holds with 3xi = 3 — x, where x is defined in formula (]68p . The limiting equation 
(I55p takes the form 

E e (x) = E (x) + / g(x, y)C(y)E e (y)dy, (76) 
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where E e {x) := lim a ^o E{x) is the effective field in the medium when the number of small 
particles tends to infinity while their size a tends to zero, and 

C(x m ) := c lm N(x m ). (77) 

Formula ([77|) defines uniquely a continuous function C(x) since the points x m are 
distributed everywhere dense in D as a — ► 0. The functions C(x) can be created at our 
will, since it is defined by the numbers c\ m and N(x m ), which are at our disposal. 

Apply the operator V 2 + k 2 to ([75]) and get 



[V 2 +IC 2 ( y x)}E e = 0, K 2 {x):=k 2 + C{x):=k 2 n 2 {x). (78) 

Therefore the (refraction) coefficient n 2 (x) is defined by the formula: 

n 2 (x) = l + k- 2 C(x). (79) 

The function K, 2 (x) is quite different from the function K 2 (x), defined in (J5J). The function 
C(x) and, therefore, n 2 {x) depend on the choice of N(x) in (|73p and on the choice of the 
coefficients c\ m in ([72]) . The c\ m is determined by j m , defined in ([57]). Since one is free to 
choose 7 m and N(x) as one wishes, one can create a desired function n 2 (x) by choosing 
suitable 7 m and N(x). 

Note that the total volume of the small particles is negligible as a — > 0. Indeed, this 
volume equals to the volume of one particle, which is -f-a 3 times the total number J\f(D) 
of the small particles in D. One has x > 0, so M{D)^fa z = 0(o~^ 3 ~^a 3 ) = 0(0") -» 
as a — > 0. 

The limiting effective field E e satisfies equations ([75]) and ([75]) . but does not, in general, 
satisfy the equation 

V • E e = 0. (80) 

Indeed, the solution to (|T6|> is unique, and it follows from (|78l) that 

V 2 ? 7 + /C 2 r ? = -V/C 2 (x)-S e , /C 2 (a;) :=fc 2 + C(x), (81) 

where 77 := V • i£ e . Since (|81|) is inhomogeneous, and V/C 2 • E e 7^ 0, in general, its solution 
7/ 7^ 0, in general. 

The unique solution to ([76]) does not satisfy, in general, the Maxwell's equation (JTJ). 
Thus, in the limit a — > one gets a vector field E e which satisfies a new equation (|76|) . 

What is the physical meaning of this new equation (j76p ? This equation is equivalent 
to equation ([78]) the solution of which is required to satisfy the radiation condition. If one 
rewrites equation ([78]) in the form: 

V xV x E e = IC 2 (x)E e + VV -E e , (82) 

then one can interpret the term VV • E e as the "current" iuifiJ in the Maxwell's equation. 
Indeed, Maxwell's equations 

V x E = iujfiH, V x H = -iu;e'(x)E + J, 
13 



imply 

VxVx£ = K 2 (x)E + iujfiJ. (83) 

Therefore the term (iu;/i) -1 VV • E e := J can be interpreted as a "non-classical" current: 
the classical current is of the form crE, where a = a%j is the conductivity tensor. In our 
interpretation tensor a is a differential operator: a = <7y = $ x < dx , so oE = VV • E. 

A different interpretation of this term may be based on the assumption that the vector 
of the induction D e in the limiting medium is related to E e non-locally: 



D e (x) = eE e (x) + \ x(x,y)E e (y)dy, 
Jd 



where xi. x iy)i the susceptibility tensorial kernel, is the kernel of the non-local operator 
defining the polarization vector in the limiting medium. Such non-local operator of po- 
larization will produce the term VV • E e if %(x, y) is a distributional kernel of the form 
x(x,y) = {u 2 n)~ l V x 8(x - y)V y . 

Indeed, Maxwell's equations in this case are: 

V x E = iuj^iH, V x H = —iujeE — iuj I x{ x i y)E{y)dy. 

Jd 

With the above choice of x( x :V)j these equations imply: 

V xV xE = JC 2 (x)E + VV-E. (84) 

This is equation (|82p. which is equivalent to (|78p . The function e in the above Maxwell's 
equations is related to /C 2 by the usual formula u> 2 /j,i = /C 2 . 

It would be of interest to understand the properties of the limiting medium from an 
experimental point of view. 

4 Spatial dispersion and negative refraction 

If n 2 (x) = n 2 (x,uj), so that uj = cj(K), then one says that spatial dispersion takes place. 
Here K is the wave vector, |K| = K. Its direction is given by the unit vector K°, which 
is oriented along the direction of the propagation of the wave front. The group velocity 
v g is defined by the formula 

v g = V K w(K), 

while the phase velocity v p h is defined by the formula 

„ =" K °, K°:=-, K:=\K\. 
P h | K | , R , ll 

Negative refraction means that the group velocity in a medium is directed opposite to the 
phase velocity. 

A sufficient condition for an isotropic material to have negative refraction was derived 
in |13j . Namely, in an isotropic medium uj = w(|K|), and 

Lun(x,uj) = cK, K = |K|, 
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where c > is the wave speed in the free space. Applying the gradient with respect to K 
to the above equation yields 

\n(x,u;)+u J ^]V K u J (\K\) = cK . (85) 

OOJ 

If Imn(x,uj) = 0, then (|85p shows that Vk^(K) is directed opposite to K° if and only if 

OTl \X LO] 

n(x, u) + u — %r-^- < 0. (86) 

ou 

To create a material with negative refraction it is sufficient to create an isotropic ma- 
terial with a refraction coefficient n(x,ui) satisfying inequality (|86p . This can be achieved 
by choosing the constants c\ m in (f68j) depending on u> so that the function C(x,u>) would 
generate by the formula 

n 2 (x,uj) = l + k~ 2 C{x,uj) (87) 

the refraction coefficient which is a real- valued function satisfying inequality (|87p . 

5 Conclusions 

A theory of electromagnetic (EM) wave scattering by many small bodies is developed. 
A rigorous reduction of the many-body scattering problem to solving linear algebraic 
system is given bypassing numerical solution of the integral equations. An equation for 
the effective field in a medium consisting of many small particles is derived. It is shown 
that by choosing large number of small inhomogeneities with special properties one can 
create a medium with a desired refraction coefficient. 
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